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Abstract 

It is showed how the Hamiltonian lattice loop representation can be cast straight- 
forwardly in the Lagrangian formalism. The procedure is general and here we present 
the simplest case: pure compact QED. This connection has been shaded by the non 
canonical character of the algebra of the fundamental loop operators. The loops 
represent tubes of electric flux and can be considered the dual objects to the Nielsen- 
Olesen strings supported by the Higgs broken phase. The lattice loop classical action 
corresponding to the Villain form is proportional to the quadratic area of the loop 
world sheets and thus it is similar to the Nambu string action. This loop action is 
used in a Monte Carlo simulation and its appealing features are discussed. 



1 Introduction 

A unified quantum theory which describes the gauge fields and the gravitation is one of the 
main goals pursued by the physicists for long time. A good candidate for accomplishing 
this comprehensive framework is the loop representation. This loop approach was intro- 
duced in the early eighties by Gambini and Trias ||ll,|^ as a Hamiltonian representation 
of gauge theories in terms of their natural physical excitations: the loops. The original 
aim of this general analytical Hamiltonian approach for gauge theories was to avoid the 
redundancy introduced by the gauge symmetry working directly in the space of physical 
states. However, soon it was realized that the loop formalism goes far beyond of a simple 
gauge invariant description. The introduction by Ashtekar of a new set of variables 
that cast general relativity in the same language as gauge theories allowed to apply loop 
techniques as a natural non-perturbative description of Einstein's theory. Furthermore, 
the loop representation appeared as the most appealing application of the loop techniques 
to this problem 0|,0. 

The Hamiltonian techniques for gauge theories have been developed during the last 
decade and they provide interesting results for several lattice models 0-|0. On the other 
hand a Lagrangian approach in terms of loops has been elusive, due mainly to the non- 
canonical character of the loop algebra. This feature forbids the possibility of performing 
a Legendre transformation as a straightforward way to obtain the Lagrangian from the 
Hamiltonian. 

In the case of non-Abelian gauge theory a major problem has been whether we can 
write a reasonably simple Lagrangian in terms of "electric vector potentials" |10|. A 



Lagrangian loop formulation will give rise to new computation techniques providing a a 
useful complement to the Hamiltonian loop studies. 

Recently, it was proposed a tentative classical action in terms of loop variables for 
the U(l) gauge theory |ll|. Shortly afterwards, we proved that the lattice version of this 
action is equivalent to Villain form for D=2+l dimensions but is slightly different for 
D=3+l dimensions |jl2|]. In fact, this action written in terms of variables directly attached 



to spatial loops seems to fail in describing all the dynamical degrees of freedom for D=4. 

This paper is organized as follows. In section 2 we show how the loops, originally 
thought up within the Hamiltonian formalism, can be introduced in a natural way in the 
lattice Lagrangian theory. We follow a different approach to that of reference [J12l: we 
show how the electric loops can be traced in the statistical lattice formulation of 4d U{1) 
theory giving rise to an expression of the partition function as a sum of integer closed 
surfaces. We interpret these surfaces as the world sheets of the electric loops. We discuss 
the connection of this classical loop action with the Nambu string action. A clear analogy is 
patent between the gauge theory in the loop representation and the bosonic and fermionic 
strings, as it was previously suspected [|l3l, [Q but never (as far as we know) demonstrated 



explicitly. The parallelism of the loop representation with the topological representation 
of the broken Higgs phase in terms of Nielsen-Olesen strings ||15[ is also pointed out. In 
section 3 we use the loop action equivalent to the Villain form for performing a Monte 
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Carlo simulation. It turns out that this action is the same as the 7 — * oo limit for the 
non-compact Abelian-Higgs theory |T^ (7 : Higgs coupling constant). 



2 The Lagrangian loop Representation 

The loop based approach of ref. 0] describes the quantum electrodynamics in terms of the 
gauge invariant holonomy (Wilson loop) 



W{^)=exp[tefAa{y)dy^], (1) 

and the conjugate electric field E°'{x) . They obey the commutation relations 

[E'^{x),W{^)] =e [six- y)dy^W{^). (2) 

J-y 

These operators act on a state space of abelian loops "0(7) that may be expressed in 
terms of the transform 

^(7) = JdM] <l\A><A\tlj >= J d^[A\i;[A]exp[-ie f Aady'']. (3) 

This loop representation has many appealing features: In first place, it allows to do away 
with the first class constraints of gauge theories. That is, the Gauss law is automatically 
satisfied. In second place, the formalism only involves gauge invariant objects. Finally, 
all the gauge invariant operators have a transparent geometrical meaning when they are 
realized in the loop space. 

When this loop representation is implemented in the lattice it offers a gauge invariant 
description of physical states in terms of kets | C >= W{C) \ >, where C labels a closed 
path in the spaim/ lattice. Eq.(^ becomes 

[Ei,W{C)] = Ni{C)W{C), (4) 

where / denotes the links of the lattice, E{1) the lattice electric field operator, W{C) = 
rizec ^(0 ^^^ Ni{C) is the number of times that the link / appears in the closed path C. 
In this loop representation, the Wilson loop acts as the loop creation operator: 

W{C')\C>=\C' ■€> . (5) 

The physical meaning of an abelian loop may be deduced from (^) and (^, in fact 

Ei\C>= Ni{C) I C >, (6) 

which implies that | C > is an eigenstate of the electric field. The corresponding eigenvalue 
is different from zero if the link / belongs to C. Thus C represents a confined line of electric 
flux. 



In order to cast the loop representation in Lagrangian form it is convenient to use the 



language of differential forms on the lattice of ref . ^7j . Besides the great simplifications to 
which this formalism lead its advantages consists in the general character of the expres- 
sions obtained. That is, most of the transformations are independent on the space-time 
dimension or on the rank of the fields. So let us sumarize the main concepts and some 
useful results of the formalism of differential forms on the lattice. 

A k-form is a function defined on the k-cells of the lattice (k=0 sites, k=l links, k=2 
plaquettes, etc.) over an abelian group which shall be R, Z, or U(l)=reals module 27r. 

Integer forms can be considered geometrical objects on the lattice. For instance, a 
1-form represents a path and the integer value on a link is the number of times that the 
path traverses this link. 

Let us introduce V is the co-border operator which maps k-forms onto (k+l)-forms. It 
is the gradient operator when acting on scalar functions (0-forms) and it is the rotational 
on vector functions (1-forms). We shall consider the scalar product of p-forms defined 
< a \ P >= J2ck tt(c)/5(c) where the sum runs over the k-cells of the lattice. Under this 
product the V operator is adjoint to the border operator d which maps k-forms onto 
(k-l)-forms and which corresponds to minus times the usual divergence operator. That is. 



< a I V/5 >=< da\ p >, (7) 

<Va\P >=< a\dp> . (8) 

The co-border V and border d operators verify 

V2 = 0, d^ = 0. (9) 
The Laplace-Beltrami operator operator is defined by 

n = Vd + dV. (10) 
It is a symmetric linear operator which commutes with V and d, and differs only by a 



^fi- 



minus sign of the current Laplacian A^ A 

From Eq.(|TU|) is easy to show the Hodge-identity: 

1 = an-iV + Vn"^a. (11) 

A useful tool to consider is the duality transformation which maps biyectively k-forms 
over (D-k)-forms. We denote by *Pco-k ^^^ dual form of the p^k form. For example, for 
D = 2, to plaquettes there correspond sites of the dual lattice, i.e. those vertices obtained 
from the original ones by a translation of vector (a/2, a/2). 

Under duality the border and co-border operators interchange: 

a = *v*. (12) 



After this digresion about differential forms on the lattice let us consider the generating 
functional for the Wilson U(l) lattice action: 



Zw = j {dei)exp{--Y,cosep), (13) 



where the subscripts / and p stand for the lattice links and plaquettes respectively. 
Fourier expanding the exp[cos^] we get 



zw= ^(rf^o^E^np(/^)e"'''^ (14) 

J —TV „ r. 



which can be written, using the language of differential forms as 



Zw=T. r (rf^0exp(Eln/„,(/5))e^<"'^^'>. (15) 

In the above expression, 6i is a real periodic 1-form, that is, a real number 9 E [— vr, tt] 
defined on each link of the lattice; V is the co-border operator; Up are integer 2-forms, 
defined at the lattice plaquettes. 

By eq. (|^) and integrating over 6i we obtain a 6{dnp). Then, 

Zw^ E exp(5:in4^(/?)), (16) 

{np;dnp=0} P 

the constraint dnp = means that the sum is restricted to closed 2-forms. Thus, the 
sum runs over collections of plaquettes constituting closed surfaces. This expression was 



obtained by Savit [0 as an intermediate step towards the dual representation. 

An alternative and more easy to handle lattice action than the Wilson form is the 
Villain form. The partition function of that form is given by 

Zv = J{d9)Y.expi-(^ II V9-2rrs \\'), (17) 

where 1 1 ... | P=< ....... >. If we use the Poisson summation formula 

/oo 
#/(0)e2-^^" 

and we integrate the continuum (p variables we get 

Zv = {2%Pvy'^^ f{d9)Y,exp{--^ <n,n>+i<n,V9 >), (18) 

where Np in the number of plaquettes of the lattice. Again, we can use the equality: 
< n, V^ >=< dn, 9 > and integrating over 9 we obtain a 5{dn). Then, 
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{n;9^=0} ^^^ 



W = {27^Pvr'''" E exp(--^<n,n>), (19) 



where n are integer 2-forms. Eq. (^) is obtained from Eq.(|T6|) in the /3 ^ oo hmit. 

If we consider the intersection of one of such surfaces with a t = constant plane we get 
a loop Ct- It is easy to show that the creation operator of this loop is just the creation 
operator of the loop representation, namely the Wilson loop operator. Repeating the steps 
from Eq.([l3) to Eq.(|l|) we get for < W{Ct) > 



< I^(a) >= |(27r/5v)-^^/2 Y: eM-^<n,n>). (20) 

n 
(dn = Ct) 

This is a sum over all closed world sheets and over all world sheets spanned on the loop 
Ct- In other words, we have arrived to an expression of the partition function of compact 
electrodynamics in terms of the world sheets of loops: the loop (Lagrangian) representation. 

There are other equivalent representations which can be obtained from the Villain form. 



First, we have the the dual representation |18| obtained essentially by using the Poisson 



identity and then performing a duality transformation. Actually, the loop representation 
for the compact U(l) gauge model is reached following this procedure but stoping before 
the duality transformation. Second, for any lattice theory with Abelian compact variables, 



the HopologicaV or BKT (for Berezinskii-Kosterlitz-Thouless) representation |19] via the 



''Banks — Kogut — Myerson' transformation [^. The BKT expression for the partition 



function of a lattice theory with Abelian compact variables is given by 

27r2 
Zvoc V exp( -<*t,A*t>), (21) 

. 9 

*t 

{d*t = o) 

i.e. a sum over closed {D — k — 2) topological forms *t attached to the cells C(£i_fc_2) of 
the dual lattice and where A represents the propagator operator. In the case of compact 
electrodynamics, *t = *m i.e. the topological objects are monopoles (particles for D=2+l 
and loops for D=3+l) and A = ^. 

Returning to the loop representation of the partition function Eq.(^) we can observe 
that the loop action is proportional to the quadratic area A2: 

Sv = -^A, = -^Y.< = -W< "'^ >' (22) 

Pv PV pg5 Pv 

i.e. the sum of the squares of the multiplicities Sp of plaquettes which constitute the loop's 
world sheet S. It is interesting to note the similarity of this action with the continuous 



Nambu action or its lattice version, the Weingarten action [^ which are proportional to 
the area swept out by the bosonic string Q 

^The relation between the surfaces of the Wilson action and those of Weingarten action has been 
analyzed by Kazakov et al in ref. |E2| . 



On the other hand, we know that in the continuum the classical action of topological 
string-like solitons, namely Nielsen-Olesen vortices, reduces to the Nambu action in the 



strong coupling limit ||T5|. The Nielsen-Olesen strings are static solutions of the Higgs 
Abelian model. A patent analogy is observed when we compare the Zoop-representation 
of lattice compact pure QED and the Si^T-representation of lattice Higgs non-compact 
QED. (We compare with the non-compact instead of the compact version because this last, 
in addition to Nielsen-Olesen strings, also has Dirac monopoles as topological solutions and 
then we have to consider open as much as closed world sheets. In ref. |Q we considered 
the Higgs compact QED which exhibits duality between both representations). 

In the case of Higgs non-compact QED model ( a non compact gauge field A^ coupled 
to a scalar field $ = l^le**^) we have t* = a*, where a* represents a 2-form which cor- 
responds to the world sheet of the topological objects namely the Nielsen-Olesen strings 
P^ and A = ^\p (M is the the mass acquired by the gauge field due to the Higgs 
mechanism). Thus, both models consist in a sum over closed surfaces which are the world 
sheets of closed electric strings (loops) and closed magnetic strings (closed Nielsen-Olesen 
vortices) respectively. The corresponding lattice actions are essentially the quadratic area 
of the world sheets in both cases. Moreover, the creation operator of both loops and N.O. 
strings is essentially the same: the Wilson loop operator ||l|] [Q. 

It is also possible to regard the connection between this two models from a different 
point of view : the Villain form of U{1) is the Higgs coupling -^ oo limit of Higgs non- 
compact QED. The standard action of Higgs non- compact QED is given by 

'5 = -fE^' + 7E^^•^(O0.+^ (23) 



It was pointed out in reference |]T6| that in the limit 7 — i> oo 

U{1) = e'^' -^ 1, 
which implies for the angular variables 

61 = 27mi, 
and so, the action (Esl) becomes 



-5- = -f E(2™p)' = -27r'/? < n,n >, (24) 

2 pes 

where Up = Vui. We note that this is just the Villain action of gauge U{\) model, Eq . (^) 
but with 27r^/3 instead of -^. 



3 Numerical Computations 



Here we shall present the results of numerical simulations carried out for the loop action 
(P^) corresponding to Villain form of U{1) model. In fact, we have simulated the dual 
action in terms of the dual integer variables *np 

S^ = -^T. *nl = -^ E(V * n,r = --1 E(E *^0^ (25) 

Pv p PV p PV p i(zp 

where *ni are integer 1-forms attached to links of the dual lattice and the integer 2-forms 
*np correspond to their lattice curl, i.e. dup = means that Up = dric where Uc are integer 
3- forms attached to the elementary cubes and according to (|12D *np = V * (ric) = V * rii. 
Note that this action implies the assignment of unbounded integer variables to the links of 
the lattice and the action is defined through the square of the ordered sum of the integers 
of an elementary square of the lattice. We have implemented a Metropolis algorithm fixing 
the acceptance ratio, as it is usual in random surfaces analysis. 

We have studied this model for different lattice sizes, imposing periodic boundary condi- 
tions. Simple thermal cycles showed the presence of a phase transition in the neighborhood 
of Pv = 0.639. To study the order of this transitions we have analyzed the energy histogram 
and we have checked for the presence of tunneling between the phases. We have not applied 
any reweighting extrapolation technique, only a direct observation of the histograms. 

In Fig. la we present the histogram of the plaquette energy density including 80.000 
measures, after discarding 40.000 thermalizing iterations, for a 12^ lattice just at the Villain 
transition point /?y = 0.639. One can observe clearly a two-peak structure. 

In Fig. lb we present the time evolution of the total internal energy during the simula- 
tion. Each point is the average over 100 consecutive iterations. This analysis shows clearly 
the presence of "tunneling" between the phases. 

Our numerical results, using this loop action equivalent to Villain form and imposing 
periodic boundary conditions {PBC), exhibit a first order phase transition signal. This is in 
agreement with the standard lattice numerical simulations of QED using the Wilson action 



(again using the standard PBC) |2^,p6l. Nevertheless, some differences between the 
simulations using Wilson or Villain actions arise. In particular, we have not observed here 
the strange persistence of the phases -and the absence of tunneling- seen in the simulations 



performed with the Wilson action. Remember that according |]2^-||2^ the first order nature 
of this phase transition seems to be a spurious effect produced by the non-trivial topology 
that periodic boundary conditions bear with. 

Typical loop configurations, obtained by intersecting the lattice with x^ = constant 
planes, are showed in Fig. 2 for two values of the coupling constant, one to the left (strong 
coupling) and one to the right (weak coupling) of the transition point. The difference 
between the two typical configurations is patent. These configurations are obtained by 
taking the last measure after thermalization with 20.000 iterations and they represent all 
the plaquettes with integer value Up ^ 0. In order to obtain this plaquettes we proceed 



as follows: first, we stored the *np 7^ of the dual plaquettes and then we performed a 
duality transformation to get the corresponding non zero Up. 

In TABLE 1 we present the spectra of plaquettes {up} (in all the lattice, not only in 
a particular cut X4 = t ) for different values of the (3v coupling. Comparing the {up} at 
/3y=0.637 and /3v'=0.641 i.e. just before and after the critical Villain coupling /?c=0.639 
one can observe an abrupt increment of non zero plaquettes. 

TABLE 1. 



Plaquette Configurations 



Pv 


n = -4 


n = —3 


71= -2 


n = -1 


n = 


n=l 


n = 2 


n = 3 


n = 4 


0.200 











15 


124386 


15 











0.500 








2 


1657 


121099 


1659 











0.555 








2 


3661 


117093 


3655 


5 








0.637 








70 


12208 


99850 


12228 


60 








0.641 








115 


15104 


93972 


15116 


109 








0.714 








299 


19573 


84639 


19639 


266 








1.000 


1 


12 


1292 


25787 


70239 


25776 


1300 


9 





2.000 


12 


569 


6695 


30141 


49474 


30337 


6630 


543 
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Incidentally, we want to remark that the apparent first order behaviour of the Higgs 
non- compact QED in the 7 — ^ 00 limit found in ref. |T^ can be clarified under this new 
perspective. In that limit it turns out a sort of "compactification" which transforms the 
original model in the pure compact QED (Villain form) we have studied here. 



4 Conclusions 

As it was mentioned, the loop space provides a common scenario for a non-local description 
of gauge theories and quantum gravity. Up to now, the loop approach was exclusively a 
Hamiltonian formalism and no lagrangian counterpart was available. A classical action for 
the Yang-Mills theory in terms of loop variables would be very valuable in its own right 
because they are the natural candidates to describe the theory in a confining phase. In 
addition, it may be useful to obtain semiclassical approximations to gauge theories or to 
general relativity in terms of Ashtekar's variables. Here we present some small steps in 
this direction which continue those ones of ref. 
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In relation with the considered case of lattice Abelian gauge theory, we observed the 
known analogy between the confining phase of lattice electrodynamics described in terms of 
electric loops and the Higgs phase described in terms of Nielsen-Olesen magnetic vortices. 

We can also ask whether the loops are no more than a useful representation or if, per- 
haps, they have a deeper physical meaning. Lattice QED exhibits a confining- deconfining 
transition, although in principle, ordinary continuum QED has only one non-confining 
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phase. However, there are studies which indicate that also there is a phase transition for 



QED in the continuum |]30[. In addition to the usual weak coupling phase, a strong cou- 



pling confining phase exists above a critical coupling Oc- This new phase could explain a 



mysterious collection of data from heavy ion collisions |^-|^. The unexpected feature 
is the observation of positron-electron resonances with narrow peak energy in the range 
of 1.4-1.8 MeV. This suggest the existence of 'electro-mesons' in a strongly coupled phase 
of QED. Moreover, a new two-phases model of continuum QED and a mass formula for 
taking account of the positronium spectrum in the strong coupling phase has been regarded 
recently |^. Thus, in principle, we can speculate about the existence in nature of abelian 



electric tubes providing a real support for abelian loops and the possibility of being on 
equal footing with the observed magnetic vortices. 
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Figure Captions. 



la Histogram of the plaquette energy density corresponding to 80.000 measures on a 12^ 
lattice just at the Villain transition point Pv = 0.639. 

lb Time evolution of the total internal energy during the simulation. Each point is the 
average over 100 consecutive iterations. 

2a Typical loop configurations, obtained by intersecting the lattice with X4 = 6, at /3y = 
0.5 (strong coupling phase). 

2b Typical loop configurations, obtained by intersecting the lattice with X4 = 6, at /9y = 1 
(weak coupling phase). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9407170vl 



This figure "figl-2.png" is available in "png" format from: 
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